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Visual Computing 



Characteristic Lines 



• Types of characteristic lines in a vector field: 

stream lines: tangential to the vector field 

path lines: trajectories of massless particles 
in the flow (non-static flow fields) 

-streak lines: trace of dye that is deposited 
into the flow at a fixed position 

time lines (time surfaces): propagation of 
a line (surface) of massless elements in time 
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Characteristic Lines 



• stream lines 

-tangential to the vector field 
-stationary vector field (arbitrary, yet fixed time f) 
-stream line is a solution to the initial value 
problem of an ordinary differential equation: 




L(0)=x 0 , ^M = v(L{u)) 

^ aU 



initial value 
(seed point x 0 ) 



ordinary differential equation 
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Characteristic Lines 



• path lines 

-trajectories of massless particles in the flow 
-vector field can be time-dependent (unsteady) 
-path line is a solution to the initial value problem 
of an ordinary differential equation: 

UO)=x 0 , dl -^ = v(L(t)j) 
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Characteristic Lines 



• streak lines 

-trace of dye that is released into the flow at a fixed 
position 

connect all particles that passed through a certain 
position (non-stationary flow) 

• time lines (time surfaces) 

propagation of a line (surface) of massless elements 
over time 

many point-like particles that are traced synchronously 
connect particles that were deposited simultaneously 
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Characteristic Lines 



• comparison of path lines, streak lines, and stream lines 




path line streak line stream line for t 3 



• path lines, streak lines, and stream lines are identical 
for stationary flows 
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Characteristic Lines 
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Arrows and Glyphs 



• visualize local features of the vector field: 

-vector itself 

-vorticity 

-external data: temperature, pressure, etc. 

• important elements of a vector: 

-direction 

-magnitude 

- not: components of a vector 

• approaches: 

- arrow plots 
-glyphs 
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Arrows and Glyphs 



• arrows visualization 



- direction of vector field 

- orientation 

- magnitude: 

■ length of arrows 

■ color coding 
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Arrows and Glyphs 



• arrows 
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Arrows and Glyphs 




• glyphs 

-can visualize more features of the vector field 
(flow field) 
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Arrows and Glyphs 



• pros and cons of glyphs and arrows: 

+ simple 
+ 3D effects 

- heavy load in the graphics subsystem 

- inherent occlusion effects 

- poor results if magnitude of velocity changes 
rapidly (use arrows of constant length and 
color code magnitude) 
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Mapping Methods Based on Particle 
.Tracina 



• basic idea: trace particles 

• characteristic lines 

• local or global methods 

• mapping approaches: 
-lines 

-surfaces 

-individual particles 
-sometimes animated 
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Mapping Methods Based on Particle 
.Tracina 



• path lines 






Mapping Methods Based on Particle 
.Tracina 



• stream balls 

-encode additional scalar value by radius 
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Mapping Methods Based on Particle 
.Tracina 




Mapping Methods Based on Particle 
.Tracina 



• stream ribbons 

- trace two close-by particles 

- keep distance constant 








Visual Computing 
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Mapping Methods Based on Particle 
.Tracina 



• stream tubes 

- specify contour, e.g. triangle 
or circle, and trace it through 
the flow 
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Mapping Methods Based on Particle 
.Tracina 



• LIC (Line Integral Convolution) 
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Numerical Integration of ODEs 



• typical example of particle tracing problem (path 
line): 

L{0)=x 0 , ^j^ = v(L(t),t) 
at 

• initial value problem for ordinary differential 
equations (ODE) 

• what kind of numerical solver? 
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Numerical Integration of ODEs 



• rewrite ODE in generic form 

• initial value problem for: 

x(t)=f(x,t) 

• most simple (naive) approach: explicit Euler method 



x(t + At) = x(t) + At f (x, t) 

• based on Taylor expansion 

x(t + At) = x(f ) + At x(f) + 0(At 2 ) 

• first order method 

• increase accuracy by smaller step size 
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Numerical Integration of ODEs 



• Problem of explicit Euler method 

- inaccurate 

- unstable 



• Example: 

x = -kx 



divergence for At > 21k 
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Numerical Integration of ODEs 



• Implicit Euler method 



Approximation of derivative 

x(t + At) - x(t) 



x(t + At) ' 



At 



Implicit timestep 

x(t + At) — x(t) + At • x(Y + At) 
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Numerical Integration of ODEs 



Midpoint method 



1 . explicit Euler step 

Ax = Atf(x,t) 

2. Evaluation of f at midpoint 



Ax At 



f m d= f x +— >f+— 



3. Complete step with value at midpoint 

x(t + At) = x(t) + At f mid 

Method of second order 
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Numerical Integration of ODEs 



• Classical Runge-Kutta of fourth order 

k l = Atf(x,t) 

k 2 = Atf( x +—,t + — 

2 l 2 2 



k 3 



= Atf 




k 4 = A/ L f(x + k 3 ,/ L + A/ 1 ) 



x(f + At) = x(t) + — + — + — + — + o(At 5 ) 
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Numerical Integration of ODEs 



• adaptive stepsize control 

- change step size according to the error 

- decrease/increase step size depending on whether 
the local error is high/low 

- higher integration speed in “simple” regions 

- good error control 

• approaches: 

- stepsize doubling 

- embedded Runge-Kutta schemes 

• further reading: 

- WH Press, SATeukolsky, WT Vetterling, BP Flannery: Numerical 



Recipes 
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Particle Tracing on Grids 









vector field given on a grid 
solve 



L(0) = x 0 



dLjt) 

dt 



= v(L (t),t) 









for the path line 
incremental integration 
discretized path of the particle 
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Particle Tracing on Grids 



• most simple case: Cartesian grid 

• basic algorithm: 

Select start point (seed point) 










r 








i 






) 








/ 


Find cell that contains start point point i 


location 



While (particle in domain) do 

Interpolate vector field at interpolation 
current position 

Integrate to new position integration 
Find new cell point location 

Draw line segment between latest 
particle positions 
Endwhile 
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Line Integral Convolution 



• a global method to visualize vector fields 
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Line Integral Convolution 



• line Integral Convolution (LIC) 

visualize dense flow fields by imaging its integral curves 

vover domain with a random texture (so called .input texture 1 , usually 

stationary white noise) 

blur (convolve) the input texture along the path lines using a specified filter 
kernel 

• look of 2D LIC images 

intensity distribution along path lines shows high correlation 
- no correlation between neighboring path lines 
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Line Integral Convolution 



• idea of Line Integral Convolution (LIC) 
-global visualization technique 
-start with random texture 
-smear out along stream lines 




Line Integral Convolution 



• algorithm for 2D LIC 

let ® 0 (?) be the path line containing the point (x 0 ,y 0 ) 
T(x,y) is the randomly generated input texture 
compute the pixel intensity as: 



I(x 0 , ;y 0 ) = J_ L k(t) ■ T(<f> 0 (0) dt 



convolution with 
kernel 



• kernel: 

- finite support [~L,L] 

- normalized 

often simple box filter 
often symmetric (isotropic) 



kernel 

k(t) 



1 J L L k(t)dt=1 
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Line Integral Convolution 



• algorithm for 2D LIC 
- convolve a random texture 
along the stream lines 
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Vector field 



path line 



Input texture 



Output image 
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Line Integral Convolution 




kernel 

Convolution k(t) 

-L 

T 
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Vector field 



f L k(t)dt=1 

L 



Final image 
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Line Integral Convolution 



• fast LIC 

• problems with LIC 

new stream line is computed at each pixel 
convolution (integral) is computed at each pixel 
-slow 

• improvement: 

compute very long stream lines 

reuse these stream lines for many different pixels 

incremental computation of the convolution integral 
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Line Integral Convolution 



• fast LIC: incremental integration 

- discretization of convolution integral 
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Line Integral Convolution 



• fast LIC: incremental integration 

- discretization of convolution integral 




Line Integral Convolution 



• fast LIC: incremental integration for constant kernel 
-stream line x 0 ,..., x„with m,n >L 

-given texture values T m ,...,T 0 ,...,T n 
-what are results of convolution: \ n -\? 

for box filter (constant kernel): 



/ = i > t 

'0 2L+1 ' i 

i=-L 



incremental integration: 



I -j- -J— 
1 j+1 j 2L+1 
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Line Integral Convolution 



• fast LIC: Algorithm 

data structure for output: Luminance/Alpha image 



■ luminance = gray-scale output 

■ alpha = number of streamline passing through that pixel 



For each pixel p in output image 
If (Alpha (p) < #min) Then 

Initialize streamline computation with X 0 = center of p 
Compute convolution /(x 0 ) 

Add result to pixel p 
For m = 1 to Limit M 



Incremental convolution for /(x m ) and /(x_ m ) 
Add results to pixels containing x m and x_ m 
End for 
End if 



End for 

Normalize all pixels according to Alpha 



a 



Visual Computing Institute | Prof. Dr. Leif Kobbelt 
Computer Graphics and Multimedia 
Geometry Processing 




Line Integral Convolution 



• oriented LIC (OLIC): 

- visualizes orientation (in addition to direction) 

- sparse texture 

- anisotropic convolution kernel 

- acceleration: integrate individual drops and compose them to final 
image 




7 anisotropic 
convolution kernel 
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Line Integral Convolution 



• oriented LIC (OLIC) 
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Line Integral Convolution 



• LIC - Line Integral Convolution 
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Line Integral Convolution 




V' a&ftj m iij 7 

Lic-Applications: length of convolution integral 



with respect to magnitude of vector field 
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Line Integral Convolution 




Lie and color coding of velocity magnitude 
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Line Integral Convolution 



• summary: 

dense representation of flow fields 
convolution along stream lines -» correlation along 
stream lines 
-for 2D and 3D flows 
-stationary flows 
-extensions: 

■ Unsteady flows 

■ Animation 

■ Texture advection 
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3D Vector Fields 



• most algorithms can be applied to 2D and 3D 
vector fields 

• main problem in 3D: effective mapping to 
graphical primitives 

• main aspects: 

-occlusion 

-amount of (visual) data 
-depth perception 
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3D Vector Fields 



• approaches to occlusion issue: 

-sparse representations 
-animation 

-color differences to distinguish separate objects 
-continuity 

• reduction of visual data: 

-sparse representations 
-clipping 

-importance of semi-transparency 
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3D Vector Fields 



• missing continuity 
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3D Vector Fields 



• color differences to identify connected structures 






3D Vector Fields 




• reduction of visual data 
-3D LIC 
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3D Vector Fields 



• improving spatial perception: 

-depth cues 

■ perspective 

■ occlusion 

■ motion parallax 

■ stereo disparity 

■ color (atmospheric, fogging) 

-halos 

-orientation of structures by shading (highlights) 
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3D Vector Fields 
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3D Vector Fields 



• illuminated streamlines [Zockler et al. 1996] 

- model: streamline is made of thin cylinders 

- problem 

■ no distinct normal vector on surface 

■ normal vector in plane perpendicular to tangent: normal space 

■ cone of reflection vectors 
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V 



T 
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3D Vector Fields 



• illuminated streamlines ( cont .) 

light vector is split in tangential and normal parts 

V R = V -(Lj. — L w ) = V -((L-T)T-(L -N)N) 

= (L • T)( V • T) - (L ■ N)( V ■ N) 

= (L-T)(V-T)-7l-(L-T) 2 Vl-(V-T) 2 
= /((L -T),(V -T)) 



Idea: Represent f{) by 2D texture 
-Access pre-computed f() during 
rendering 



normal spac 
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